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1 Introduction 



The presence of dark matter in the universe has been invoked to explain 
a number of phenomena, the rotation curves of galaxiesHj-lQ the apparent 
mass of galactic clusters [|, |[ and the observed expansion of the universe. 
Many attempts to detect this conjectured matter are under way but at the 
moment there is no direct observation of it. In this paper we present a model 
with a modified version of gravity that can accommodate these points in a 
relatively straightforward way. It does so without the necessity of postulating 
dark matter. The theory, which is geometrical in character, is perfectly 
covariant, locally Lorentz invariant and satisfies the Equivalence Principle. 
It also satisfies all the tests of General Relativity on the scales at which they 
have been carried out. 

The theory is a bi-metric theory. Such theories have a long history [|^, p, |^ 
and have recently been used as a way of realising Variable Speed of Light 
(VSL) theories [|10]-|Q. An earlier incomplete version of the theory presented 



here was motivated in this way Bi-metric theories ran into difficulties 
when, as then formulated, they appeared to be inconsistent with the tests 
of General Relativity However the theory proposed here does meet 



all those tests E^, |2B|, ESI, E3, ER EBI. Of course VSL is built into bi-metric 



theories. This will be important in applying the model to early universe 
studies. However in this paper we concentrate on those implications of our 
theory for the gravitational interaction of matter that provide an explanation 
of the current state of the universe without recourse to dark matter. 

Our proposal is to introduce into the space-time manifold, two vierbein 
bundles. Each bundle supports its own metric. One is associated with matter 
and the other with (underlying) gravity. The matter vierbein can be strained 
and scaled relative to the gravitational vierbein. The dynamics of the theory 
includes this straining and scaling as dynamical degrees of freedom. The 
justification for introducing these new effects is ultimately in the results that 
emerge. However if dark matter is not present some modification of gravity or 
dynamics is essential. We choose to modify gravity in a way that permits the 
introduction of a galactic scale, something that is impossible with standard 
General Relativity. 

The equations of motion are obtained from an action comprising three 
parts, a gravitational term Iq-, of the standard curvature form, a matter term 
Imi of the standard form based on the matter metric, and a linking action 
II, that depends on the variables that determine the relationship between 
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the two vierbein bundles. The full action, 7, is the sum of all three terms, 

I = Ig + Il + Im . (1) 

Each term has its own gravitational coupling constant with the dimensions 
of Newton's constant, Gn ■ In developing the theory we show how these 
different constants are related to one another and to Newton's constant. 



2 General Structure 

The geometrical character of the theory is clearly revealed by its formula- 
tion in the vierbein formalism. The local Lorentz invariance in both the 

gravitational and matter vierbein frame is explicit throughout as a gauge 
invariance. We introduce a vierbein bundle appropriate to gravity, {e^a}, 
with the associated metric 

where raising and lowering of a-indiccs is carried out with the standard 
Lorentz metric riah — {1, —1, —1, —1} ■ The inverse vierbein is {e"'^} so 
that 

e^e'^^^Sl , e'^^e,,^5l , (3) 

and 

g^'^^e'^'e,^ . (4) 

The vierbein associated with matter is {e^a} a-nd the raising and lowering 
of a-indiccs is by means of the Lorentz metric, 77^5 = {1,-1,— 1,-1} . The 
associated metric is 

The two vierbein bundles are related by a local linear transformation 

e^-a^e^aM'^-.e^ , (6) 

where the matrix M is an element of >S'L(4, R) and the local scaling is intro- 
duced through the factor . The determinants associated with the volume 
elements of the bundles are J and J where 

J = det{e^"} , and J = det{e^"} . (7) 
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We have therefore J = Je^"^ . We denote the inverse matrix by ^ so that 

M%M% = 5» , M\M% = 5f . (8) 

The vierbein connections and associated coordinate connections are de- 
fined so that the vierbeins are covariantly constant in the appropriate way. 

D^eua = df,eya + ^(^^b " T^^CAa = , (9) 

and 

D^e^-a = d^e^a + t^^s ^^ub - Tj^CAa = . (10) 
The requirement that rjab and rj^i be covariantly constant imphes that uJ^ab = 

—u^ba and u^^h = ~^ijha ■ 

It is convenient to define a covariant derivative of M that includes both 
the right and left vierbein connections, 

i^^M^^^^M^+^x.-^^'^;. ■ (11) 

However further differentiation requires the coordinate connection, F^^^ . The 
covariant deriviative Dfj^ can be defined and extended in a similar way Its 
effect on Af is the same as that of but a second differentiation must use 
the coordinate connection . 

Gravitational curvature tensors appropriate to each of the bundles are 
defined so that 

and 

V-a^R-ab^y , (13) 

where includes the left vierbein connection field, uj^ab but not the coor- 
dinate connection, F^^, . Similarly includes the only the right vierbein 
connection, uj^i . We have then 

(14) 

with a similar definition for RaSf^u ■ It follows that 

[D„ D^] M'^, = R'^ ,^,M' , - M» -,R' - 2Ci,D,M'^ , , (15) 

where 

ci. = ^ (rj. - J . (16) 

The quantity C^^ is the torsion tensor in the coordinate basis. It is not in 
general zero. 



K — RabiivV^ 



(12) 
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3 Gravitational Action 

The gravitational action has the standard form 



G 



J:^ji^.JR , (17) 



where 

= e^V'^i?^^^, , (18) 

and G is a couphng with the dimensions of Newton's constant, Gn ■ The 
vierbeins and the connections are treated as independent variables. If we 
vary the former then 

<^/g = ^ / d'xJ5e,, [e'^^'R^ ^ - U^^r) . (19) 
The variation of the vierbein connection yields 

SlG = -Y^Jd'xJe''>^e"'{Df;Suj,ab-DiSuj^ab) • (20) 
Switching to the full covariant derivative we get 

SIg = / d'^xJe'^'^e'" {D^Su^ab - DJu^ab + '^C^.^uxab) ■ (21) 

Using the result 

JD^V = d^iJV) + 2JC^^Vi' . (22) 
we can integrate by parts and obtain finally 

SIg = / d'xJe^^e'^ {cldl - 0^6; + C;^) Sco^^b ■ (23) 

If there is no other interaction in the theory we can deduce from the vanishing 
of these variations that C^^ = 0, and 

R"^ ,, — —5ZR = , 
the standard equations for matter less gravity. 
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4 Linking Action 



Proposals for constructing an action for the linear transformation relating 
the bundles have involved parametrizing it in terms of a vector or scalar field 



16| , \17\ , pig , p!^ , pO| . Our proposal treats all the degrees of freedom inherent 
in the transformation. This is crucial for the structure of the theory. 

Because scaling commutes with the other elements of the group of linear 
transformations, it can be treated separately. The Lagrangian for can be 
chosen to be proportional to g'^^df^cpdiycj) . The remaining degrees of free- 
dom are represented by the matrix M which lies in the non-linear manifold 
SL{4, R) . A natural way of constructung a Lagrangian for such a theory is 
to invoke the mechanism of the non-linear sigma model and express it as a 
quadratic form in the derivatives {D^M)M~^ with the appropriate structure. 
The use of the covariant derivative guarantees the local Lorentz gauge invari- 
ance relative to both bundles and the presence of guarantees that the 
derivative is a proper element of the tangent space to the manifold 5'L(4, R) . 
We take as our action 

II = J d'xJg'^'^Trij^j,) + J d'xJg^''{d,<j>d,<P) , (24) 

where F and F' are new gravitational constants with the same physical 
dimensions as G . The matrix valued current is given by 

3, = {D^M)M-^ , (25) 

or more explicitly 

3f = {D,M'^-^M'' . (26) 

It is also convenient to define an alternative version of the current, ap- 
propriate to the barred vierbein bundle, as 

jf = [m-'D^mY' = M\D^M-' . (27) 
We also include in "mass" terms of the form 

1 r .A , 1 



(28) 
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For simplicity we have chosen the mass parameter m to be the same in both 
these these additional terms. For the choice 7 = 8, the action vanishes when 
M represents a Lorentz transformation. Departures from this value introduce 
a cosmological constant term in the action. By construction these additional 
terms clearly maintain local Lorentz invariance in both vierbein bundles. The 
mass terms are crucial for the effectiveness of the theory because we identify 
the galactic distance scale with . For effects on a scale much less than 
this therefore m will be viewed as a small parameter that can be neglected 
in certain circumstances. 

We treat the vierbein, which enters through g'^", the matrix, M, and the 
connections cUi^ab ^-nd a)^^^ as independent variables. The result for the linking 
action from the vierbein variation is. 



SttF 



SttF 



SttF' 

From the left vierbein connection we have 



^— j d^xJ5e,,e-^02 _ (29) 



= ^ / d^xJdu^ai.f^'' , (30) 
and from the right vierbein connection 



5It 



^ J d'xJdu;^,if'-^ . (31) 



SttF 

On varying the matrix M we obtain 

i2 



6Il = J d'^xJ{6MM 



-l\ab 



4 

(32) 

where the square brackets [■ ■ ■] indicate the traceless version of the quantity 
contained within them. The quantity 5MM~^, being an arbitrary element 
of the SL{4:, R) Lie algebra, is sufficiently general to identify the other factor 
in the integrand. Finally, on varying 0, we have 

= / d'xJS<P {g^^ (D,d^4> - 2Cld,(t>) + m^(t>) ■ (33) 
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5 Matter Action 

We assume that matter is propagated in the vierbein background {e^s} ■ 
This seems a consistent approach since it imphes that matter behaves in 
a conventional way in relation to the gravitational field it experiences. In 
particular the Equivalence Principle is satisfied. However the theory docs 
change the relationship of this observed gravitational field to the distribution 
of matter density. We have 

SIm d^xJSg^.f^'^ , (34) 

where T'^'^ — T"^ is the symmetric energy momentum tensor for matter. 
Since 

it follows that 

5Im = - j d^xJ5e^-ae, ^f'^'^ . (36) 

However 

5e^^ = 5e^M- + e^„5M<^ ,e<^ + 50e^, , (37) 
so the variation of the matter action takes the form 

5Ij^^- J d^xJe^^ {Se^cT"" + Tr(5MM-^U) + (50T) . (38) 
where f = Qf^^.f'"', 

U^" ^e^^M^ -^e^'^f"" -^r]°'^f , (39) 

and 

T"^ ^ M^e^^f"" , (40) 

T\ = e^.T'^' = e,se, '^f'^" = f'^^gx, . (41) 

If we adopt the convention that barred quantities, that is those appropriate 
to the gravitational background of the matter, have spatial indices raised and 
lowered with the barred metric we can define 

f\ = f''''gx. . (42) 
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Hence we get the simple seeming result 

T\ = T\ . (43) 

However it is important to recall that 

T"^ = T^^g^^ ^ f"^ . (44) 

In fact T'^^ is not necessecarily symmetric. This does not cause any difficulty. 
Note that T = T^" = f . 

6 Equations of Motion 

We obtain the equations of motion by requiring that the variation of the 
total action is zero. The result is 



SttG V " 2 f J SttF V ' 2 ^ 

1 TD^ 

5" (M" ,M « + M\M^ - 7) 

i2 



8nF 



1 / 1 \ 1 rp'^ 

^[9''''d^<pd,<P)--6;{g'^d,<pd^<f>)y^—6l^^ = e^^f%(45) 



Eq (|45|) is a generalization of the standard equation of General Relativity. The 
new straining degrees of freedom are controlled by the following equation, 

^ [D,ji - 2Cl3l + ^ [M, -^M,, - + e'^'U"^ = . (46) 

The generalization of the "no torsion" rule in General Relativity is repre- 
sented by the next two equations, 

r^''"^ = ^e^'e'^ {Cl5: - ClS; + 0%) . (47) 

= . (48) 
Finally the local expansion is controlled by 

^ {D,d,<j> - 2Cld,<j>) + m'<j>) + e"^f = . (49) 
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Eg (^31) implies that 



StiG ^ 8nF SttR 
1 



8nF 8 



6l{M\M,-^ + M\M,'^-^) 

+ ^!!lsi(j)^ = e'^(f\^--6:f] . (50) 



Although rather complex these equations are surprisingly susceptible of anal- 
ysis as we show below. 



7 Bianchi Identity 

Just as in standard General Relativity it is necessary to check that the theory 
satisfies the integrability conditions associated with the Bianchi identity. In 
the presence of torsion this is changed to the following 

where ; /i indicates the covariant derivative . In the contracted version it 
becomes 

{r' . ~ \^'.R) ^ = R'"" p.c;. + 2i?^' pCP, . (52) 

If we take the covariant divergence of the left side of eq (^5D and make use of 
the equations of motion together with the identity 

J? = e'^^ (f - e/ - r/'^^S,^ , (53) 

we obtain the result 

D^f\ = 2C:xf-2Clxf\ . (54) 

It is easily checked that this is equivalent to the standard conservation law 

D(5)f% = , (55) 

where D^^^ is the covariant derivative formed from the metric connection 
arising from g^y . 
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8 Connection Structure 



In standard General Relativity the assumption that matter couples to grav- 
ity only through the metric means that torsion plays no role in the theory. 
In the theory presented here we make the same assumption about matter. 
However because of the extra complexity of the theoretical structure, torsion 
is in general not zero. Nevertheless the equations of motion do permit the 
connection structure to be elucidated in a straightforward way. 

The starting point of the analysis is eq(^) which reveals that the an- 
tisymmetric part of j^ab vanishes. It allows us to express cu^^S terms of 
the other dynamical variables and hence eliminate it from the equations of 
motion. 

We now express jf^ab in terms of the symmetric part of j^^ab as follows 

J,,, = %,-,M^ , = M„ %{,,5}M^ , . (56) 
We can write this more explicitly in the form 

(57) 

which shows that j^ab is a linear function of cj^af, . 

The torsion tensor is also a linear function of uj^ab ■ The covariant con- 
stancy of e^a implies that 

ri. = e'^' (d^eya + oo^a'eub) 

Hence ^ 

C^u = 2^"^ {df.e,a + uj^, *e,fe - d^e^a - uj,, ^e^^) . (59) 

If we define 

Cx^lu = gxaC^y , (60) 

and 

= "e^ ''u^ab , (61) 
then we obtain uj^ab as a linear function of the torsion, 

f^/iAi/ = C\nu — CufiX + CfiXu + ^fiXu , (62) 



(58) 
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where cD^ai/ is the metric version of uo^x^^ and is given by 

i^iiXu = ^ (e^ "df.exa - "d^e^a + d^gx^ - dxgu^i) ■ (63) 
We note that eg (^7|) relates the torsion tensor C^^ hnearly to j^ab ■ 



Cafii^ + Cxudua — CxyQ^a — ^^afw ; (64) 



where 
and 

It follows that 



G 

Xa^^u = j;Ja[t,,u] , (65) 



Ja,,u = "e^ ^jaab ■ (66) 

Caiiu = —X„^y — -Xx^Qy^j + -XxyQ^a ■ (67) 

Finally we have the linear relation 

^pl\u = —Xxfiu + Xy^j^x — X^xu — Xlxdnv + Xl^g^x + ^^Xv , (68) 

that determines cj^ab in terms of the other dynamical variables. 

For future reference and to show that the above equation can take a 
simple form in special cases we compute jxab when the matrix M takes a 
diagonal form, namely 

M% = AA" and M\ = A'^^ . (69) 

We enforce an obvious correspondence between the values of the symbols a 
and a to give meaning to the 6 symbols. We can compute j^ab from eq(| 
to obtain 

JXab = ^Vab + lcOxab(l-§] , (70) 



with the simple result 



^AM = Ju;Aa. (2 - ^ - ^) . (71) 



Used in conjunction with eq(|68D this leads to an easy evaluation of the vier- 
bein connection. 
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9 Weak Field Limit 

In order to apply it to planetary motion, galaxies and galactic clusters it is 
appropriate to examine the theory in the limit of weak gravitational fields. 
Since we are considering relatively local objects in an effectively fiat back- 
ground wc will choose the parameter 7 = 8 in order to set the cosmological 
constant to zero. 

The weak field limit has the form 



where 

Similarly we can set 
where 



C/ia — + h/jM , (72) 



^^°ua^*"°l ~ ) (73) 



g _ ^ g(0) , h _ (7A) 

e%e^%~^^V,. . (75) 



The connections uj^ab-i '^^al ^^'^ ^^e scaling field are first order quantities. 
The matrix M has the form 

M% = M(°)% + m% , (76) 

where M(°) is a constant Lorentz transformation. We have then 

^(0) -a ^ ^ ]^{0)^ -a _ ^77^) 

Wc can use M^^^ and its inverse and e^*^^ and e*^°^ to convert superfixes and 
suffixes between the various bases. For example we have 

m\ = m'' ^M^°^ " ft . (78) 

The requirement that det M — 1 implies that 

m\ = m% = m% = . (79) 

The relationship between e and e implies that 
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or 

h^i, = h^y + m^j, + </)?7^j, , (81) 

together with corresponding equations in other bases. 
To lowest order 

Rabiiv = d^LOyab — d^UJ^ab ■ (§2) 

Hence 

= e^'^'^eW^^K^M. , (83) 

so that 

= 5.^. - 5.^, '^^ , (84) 

and 

R = 2d^u,^"' . (85) 
Again in the lowest order approximation 

jfiab = df.niab + UJ^ab " i^fiab ■ (86) 

If we convert to the coordinate basis we have 

j^,XT = d^mxr + ujf,XT - ^^^lXt ■ (87) 

We can also evaluate . In this lowest approximation it coincides with . 
From eq(|4^) we have 

JM[A,r] = JM[A,r] = . (88) 

That is 

dam[x,T] + t^aAr - ^^ctXt = . (89) 

From eq(p7D we see that the torsion in the gravitational vierbein bundle 
vanishes. Explicitly we have 

C',. = l{d,K'-d.h/ + uj^\-ujy\)=0 . (90) 

EqdD yields 

r]f"'df, {d^mxr + uj^Xt - u^aXr) + m'^m{x,T} = StcFUxt , (91) 

where 

Uxr = Txr-^VXrf , (92) 
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and we have assumed that T^,^ and hence U^i, is a first order quantity. Making 
use of eq(|89D we obtain the result 

V'^'d^d^m^x^r} + m'm^x^r} = SnF {f^r - ■ (93) 



Because we set 7 = 8 in eq(E5|) we obtain 



RaX - ^VaxR = SnGf^x . (94) 

In the present approximation m[^^j,] can be removed by gauge transfor- 
mations of the form 

l^aAr i^aXr + f^a^Ar , aud U^Xr ^ ^^aXr + f^a^Ar • (95) 

We can assume therefore that in this approximation rn^^^^-^ vanishes. There- 
fore m^,^ may be assumed symmetric. It satisfies 

{d^ +m^)m^, = -87rF{f^,-^ri^,f) . (96) 

The gauge invariance referred to above means also the we are free to choose 
h^i, to be symmetric with the result that /i^,^ is also symmetric. Under these 
circumstances we can solve eq (|90|) to yield 



i^uXf, = d^Kx - dxKf, . (97) 

Eq (|45|) now implies that 

(98) 



We now refine our coordinate system by choosing the harmonic gauge. 

g^^Tl = Q . (99) 
In the lowest approximation it yields 

1 

— t 

2 

The equation of motion then becomes 



d,h\ = -dxh\ . (100) 



( h^^a - ^Vf^ah^ r ) = -SttGT^, , (101) 
15 



or 



a V = -s^G- (7;, - ^77^,r) . (102) 

The equation for (j) is 

(a^ + m^)^ = -87rF'f . (103) 

To demonstrate how the theory matches up to the standard tests of Gen- 
eral Relativity we consider its implications for phenomena on a scale such 
as the solar system that is much smaller than the galactic scale. For such 
applications we can set m = . We will return to the problem with the 
galactic scale parameter later. 

By combining the above equations we find for the massless case 

= -8tt{G + F)f^, + Stt + ^F- F') rj^^f . (104) 
By making use of the result da-T^ — 0, we can show also that 



where 



d^h\ - ^d,h^ ^ = , (105) 



+ — F' 

2 4 . (106) 



a - AF' 

Assuming then that the relevant fields are sourced by the energy momentum 
tensor we can maintain the condition 

dJi\-idJi\ = Q . (107) 

For arbitrary F' this gauge condition on h^i, is different from that on h^^ . 
However the choice 

F'^-\f , (108) 
is of special interest. We have then ^ = | . As a result h^t, satisfies the same 



harmonic condition as h^i, . 



dji^ u - ^ = . (109) 



1 
2 

The wave equation for /i^,^ also takes a significant form and becomes 

1 



a2V = -87r(G' + F)(7;,--77^,f) . (110) 
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It follows that for this choice of F' weak field gravity is related to the matter 
distribution exactly as in General Relativity if we set Newton's constant 
Gn = G + F . Hereafter we will assume that F' has this special value and 
that Newton's constant is obtained in this way from the theory. 

For example in a static situation where d'^ = — V^, Too = T = p (the 
density of matter) we have 

V'/ioo = 47rG,vP . (Ill) 

Of course we can interpret Jiqo as the gravitational potential experienced by 
a material particle. In addition the spatial part of the metric satisfies 

V%j = A7rGMp5ij . (112) 

This is precisely the form for the spatial metric to yield Einstein's prediction 
for the deflection of light and to satisfy the time-of-flight measurements of 



radio signals |18, 19 



The remaining solar system scale test is the precession of the orbit of 
Mercury. This requires a higher order correction than the Newtonian ap- 
proximation of the weak field limit. For reasons of space we do not present 
this calculation here but we have checked the consistency of the theory on 
this point by examining the asymptotic behaviour of the Schwarzschild-like 
solution. The result is that the refinement of the Newtonian potential that 
produces the precession is correctly given by the theory. 

It is clear from the wave equation for h^iy, eq (|110|) , and the associated 
gauge condition, eq( |109D , that the gravitational waves emitted by a time de- 
pendent matter distribution will be exactly the same as predicted by General 
Relativity. The detection of these waves by ordinary matter will also be en- 
tirely conventional. It is reasonable to conclude that the observations of the 
slowing of a binary quasar and its conformity with the predictions of GR will 
be reproduced in our theory |^ ^ p6| . 

There are circumstances in which the bi-metric theory could show differ- 
ences with General Relativity. These would occur were there to be a form of 
matter that coupled directly to the metric g^^ . Such matter would act as a 
source for and react to the unbarred metric. In the weak field limit it would 
therefore behave as if Gn = G and if mixed with other matter would imply 
that the Equivalence Principle did not hold. Part of our theory of ordinary 
matter and gravity is that such anomalous matter is not present. This is the 
assumption we use throughout this paper. 
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10 "Dark Matter" and the Galactic Scale 

In order to apply our theory to objects of galactic or extra-galactic size we 
restore the "mass" parameter m . Our hypothesis is that is a length of 
the order of 30 kpc. Assigning such a value to achieves the purpose that 
the theory exhibits three regions of length scale, namely (i) 0-1 kpc in which 
conventional Newtonian gravity holds sway and for which the theory of the 
previous section is relevant, (ii) 1-100 kpc, a transition region appropriate to 
galactic dynamics and (iii) above 100 kpc in which Newtonian gravity with 
an enhanced coupling appears. 

To analyse the gravitational effect of the galactic scale we set 

h^^ = + h'^^ , (113) 

where 

= + (^r/^^ , (114) 
while h^y satisfies eq (|102| ) h'^^ satisfies 

(a2 + m2)/i;^ = -87rF(f^,-ir7^,f) . (115) 

Therefore while a highly localised matter distribution of mass M, at the origin 
yields 

/loo = -GM- (116) 
r 

it gives rise to 

-mr 

K, = -FM—- , (117) 
for the remaining part of the metric. We have then 

S„ ^ (l . |e-™) . ,118) 

It follows that the effective gravitational coupling for mr « 1 is Gat = G+F 
while at distances beyond the galactic scale for which mr >> 1, the effective 
coupling is G . Our hypothesis is that Gat < G so that gravity is weaker at 
short distances than at long distances. This is the basis of our explanation 
for "dark matter" . For this to be true we clearly must have F < . We set 

F = -eG , (119) 
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with e > 0, so the above gravitational potential becomes 



(1 — e) r 

More generally for a matter distribution p(r) we have 

= - / AV(.0 (f_"^/^-;i;!"''' . (121) 

We believe that this choice of a negative sign for F is in fact the "natural" 



choice and stabilises the theory against tumbling of coordinate frames p2| . 
We will pursue this analysis in a future publication. Here we explore the 
more immediate physical implications of our assumptions. 

10.1 Galactic Clusters 

For objects much greater than galactic size such as galactic clusters, there 
is ample evidence that the observed matter is insufficient to account for the 
gravitational potential inferred from applications of the virial theorem to the 
motion of galaxies in the cluster || . Currently the picture of such clusters 
is that 5-10% of the mass is galactic in origin with a further contribution 
from hot gas. The bulk of the gravitational potential is accounted for by 
dark matter. It is also significant that a detailed analysis of clusters suggests 
that the dark matter distribution follows that of the visible matter ||^. The 
clusters appear to be condensed versions of the local background. 

The explanation of this effect is straightforward in our model. The ap- 
pearance of dark matter is simply the consequence of an enhanced gravita- 
tional coupling of visible matter. The parameter e represents the dark frac- 
tion of apparent matter. On the basis of the above observations we should 
expect e ~ 0.9 — 0.95 . Note that because our extragalactic dynamics is still 
Newtonian, although with an enhanced coupling, it is still possible to apply 
analyses of galactic clusters that rely on the virial theorem for an inverse 
square law of force between galaxies 0. 

The inverse square law is also crucial for the argument relating the motion 



of the Local Group to optical flux due to galaxies The close alignment 
of the motion of the Local Group through the CMB with the net optical 
flux confirms that fiuctuations in the distribution of visible and dark matter 
are closely correlated. Disparities in the visible and dark matter distributions 
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would tend to destroy this alignment. In our theory of course, the distribution 
of dark matter is identical to that of visible matter of which it is merely a 
reflection. There is no room for any bias between visible and dark matter 
fluctuations. 

We will see later that the re-interpretation of dark matter as an enhance- 
ment of the Newtonian constant at large scales extends also to the dynamics 
of the expanding universe. 

10.2 Galactic Rotation Curves 

On a galactic scale we are concerned with the rotation curves of galaxies. 
The apparent asymptotic flatness of many of these curves is usually taken as 
the most direct evidence for dark matter Our explanation rests on the 

modification of the gravitational interaction described above. It turns out 
that by choosing a value for e in the range suggested by large scale dynamics 
it is possible, using a disk model, appropriate for certain spiral galaxies, to 
compute rotation curves that exhibit the features of observed rotation curves. 

As an example we construct a model galaxy with a thin disk of surface 
density 

a{r) = ade-""' . (122) 

The mass of the disk is = 27Tad/a'^ ■ The gravitational potential in the 
plane of the galaxy is 

^(^) = / dh'e-'^'^'-'-^ (l - ee-'^l'-'-'l) , (123) 

1 — ej |r — r'l^ ^ 

and the tangential rotational velocity, v, is given by 

= r. V^(r) = ^^^^ / dh'e-'^^'-'r.r'-^ (l - ee-"'^^-^'^) . (124) 
1 — e J |r — r'l^ ^ 

The galaxy NGC 3198 has a well measured rotation curve and it is ac- 
cepted that a simple exponential disk model gives a good account of its 
luminosity distribution [Q. We assume a mass distribution with the same ex- 
ponential shape. In Fig |I] we show the resulting rotation curve where we have 
chosen the galaxy parameters to be = 0.38kpc~^ and = 2.9 x lO^^M©, 
these are quite close to the values of a previous analysis @]. The theory 
parameters are chosen to be e = 0.937 and m = 0.035 kpc~^ or m^^ = 28.6 
kpc. We stress that these parameters are not a best fit but merely the result 
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of eyeball exploration. However it is not easy to reduce the value of e value 
by much and obtain a convincing shape. It is encouraging that e does lie in 
the range we anticipated from our discussion of extragalactic structure. 

The Fig Q also shows the standard Newtonian curve for the assumed mass 
distribution. The difference between the two curves is normally attributed 
to dark matter. Here we achieve the same effect by means of our modified 
gravity theory. Fig |^ shows the extrapolation of the curve to larger radius. 
There is clearly a beginning of a fall-off around 90-100 kpc although the 
there is still substantial rotational velocity out to 400 kpc. The rotation 
curve becomes Newtonian but with an enhanced gravitational coupling 



^"V^ • ^^^^^ 

Of course the challenge to our theory is to fit all galactic rotation curves si- 
multaneously with common values for m and e and plausible masses and and 
mass distributions for the galaxies. This is no easy task since many galaxies 
have a more complicated structure than NGC3198 and detailed modelling 
will be required to determine the adequacy of our theoretical predicitions. 
We intend to pursue this task in the future. However one prediction of the 
theory does not require such a detailed attack. We predict that all galactic 
rotation curves will fall away in the range 100-200 kpc from the galactic cen- 
tre. Measurements in this range and beyond would be a direct test of our 
ideas. 

Our proposal of a galactic scale determining the shapes of rotation curves 
is in conflict with the construction of a proposed universal scaling curve 



for galaxies [^. This approach takes the optical radius of the galaxy as 
the significant length scale and uses a statistical approach to establish well 
defined rotation curves. We are proposing a model in which the ratio of the 
optical radius to the galactic scale is the significant quantity for establishing 



the shape of rotation curves. The statistical approach deployed in ||28|, which 
superimposes results from galaxies of different sizes, is therefore not open to 
us. It may well be however that the detailed modelling required to test our 
theory will provide an equally good description of actual galactic rotation 
curves. More importantly, if a common scale for the underlying structure 
of all galaxies were established it holds out the prospect of better distance 
measurements and hence the possibility of measuring the Hubble constant 
with more certainty. 
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10.3 Expanding Universe 

Finally we wish to show that when the matter density is small as in the 
present epoch, the expansion of the universe is controlled by the standard 
equations of General Relativity but with Newton's constant, Gn, replaced 
by the enhanced constant G . In other words the dark matter component 
required to relate the visible matter to the observed Hubble constant is again 
supplied by the enhancement mechanism. Of course a complete reconcilia- 
tion of the equations with observation requires a contribution to the energy 
density from the Cosmological Constant. This necessity can be catered for 
also in our theory. 

To apply our theory we take the usual starting point that the spatial sec- 
tions of the universe are isotropic and homogeneous. In the present approach 
we meet these requirements by assuming that there exist basis vierbein fields, 
E^a such that EtaEtbt]"'^ — 1 and EtaE^hi]"'^ = for 7^ t . We choose our 
coordinates so that Eta is constant in space-time and E^a depends only on 
coordinates on the spatial section when \i ^ t . We finally fix the nature of 
the space-time model by choosing structure constants /^ca so that 

d^E,,-d,E^, = -U^E^''E^'^ . (126) 

Of course /bca vanishes if any of the suffixes is (timelike). For a spatially fiat 
universe all structure constants have the value zero. For a universe of positive 
spatial curvature the spatial constants /jfcj = 2ejfcj, where 2 is a convenient 
normalization, and for a negatively curved universe /j^i = bijrik — Sikfij where 
Uk is an arbitrary unit three- vector. 

We construct the expanding universe by choosing the gravitational vier- 
bein to have the form 

ei,a = AaE^ , (127) 

where 

Ao = l and A,, = A{t) , (128) 

A similar structure is maintained for the matter vierbein by requiring the 
transformation matrix M to be diagonal. That is (we equivalence the labels 
a and a in the obvious way), 

ei,a = e'^Aaei,a = AaE^a and A„ = A^A^e^ , (129) 

where Aq = A(t) and Aj = As{t) ■ The requirement that detM = 1 implies 
that A5 = A-3 . 



22 



Using the above information and eq(^) we find that 

'^iOj = —l^ijO = ' (130) 

and ^ 

i^ijk = ^■^'^'^ ~^ ■^''^^ ■^''^''^ ' (131) 

All other components vanish. A similar structure is found for Uabc ■ From 
eq(|7I|) we obtain 

where 

G ( A2 K\\ ^ ^ 

and 

^ijk = (^ijk ■ (134) 

From eq(|14D we find for the relevant components of the curvature tensor 

Roo = Adt(i] , (135) 



and 



A ' \E 



1 „ /i\ / i ' 



~4A^ (/fcji + fijk + //fcj) (/ijfc + fkii + fkii) 
If we evaluate this expression for the three curvature cases we find 

where conventionally = 0, ±1 according as the spatial curvature is zero, 
positive or negative. We also have from eq(|70|) 
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and 



JiabJ j 



ba 



1 / A 



2 \AE 



(139) 



If we denote the matter energy density and pressure by p and p respec- 
tively, then from eq(|^) we obtain 



A 



AG /A' 



G 



+ 



A '\E 
Gm^ 



3F \ A 



F 8 



(140) 
AnGe^'f' {p + 3p) , 



and 



A 

Gm'^ 
'F~8 



E 



2E 



A 

AE 



2k 



(a^ + A-^ + 3A| + 3Ac 
From eg (|46|) we find 
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G 



(141) 
AnGe^^ip-p) . 



dt 



'a: 



+ 3 
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A 



(142) 
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(A^-A 



-2 



A| + A- 



-67rFe^'^ (p + p) . 



Finally from eq(|49|) we obtain 



3^ i 2 
+ —0 + m 



TiF'e^'^ {p-3p) . 



(143) 



Eq (|141|) and eq( |142| ) can be combined to eliminate the second derivative 
in t . 
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AE 
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2G /A' 



A^ 9F\A 
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6F'^ 



(144) 



A^(A^ + A- + 3A| + 3A--7 
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If we multiply this equation by and differentiate with repect to t we can 
use the above equations of motion to deduce that 

P= (^-3^-30] (p + p) . (145) 



This is easily re-expressed as 

"a / ' ^^M a 



or 



dt(pA^)+pdt{A^) , (147) 

where A = e'^AA~^ is the cosmic radius parameter appropriate to matter. 
Eq( |147|) is therefore the standard equation for the conservation of the energy- 



momentum tensor in this special case, as we should have expected. 

The above equations are rather complicated but can be simplified if we ask 
how they might be applied to our expanding universe in which the pressure 
vanishes and the matter density is low and getting lower. In the absence of 
matter eq( p.45| ) has a solution for which = 0. For weak density we assume 



that there is a solution for which ~ 0(p) . We have then as a leading 
approximation 

3A 

+ —0 + = -SnF'p . (148) 

The form of this equation suggests that, ignoring transients, the solution is 
to a good approximation 

= -p = -p . (149) 

We expect then that will remain 0{p) . If we set A = and omit all terms 
0(^^) then ~ 1 and eq( |143| ) becomes 




e + 3^e+ |"^'-8(^) le = -6vrFp. (150) 

Again we expect there to be a solution ^ = 0{p) . If now we neglect all terms 
O(p^) we obtain the equations 

= ^{p + Pcc) , (151) 
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where the energy density associated with the cosmological constant, pcc is 
given by 

and 

. (153) 



These are the standard equations for the expanding universe but with New- 
ton's constant Gn replaced by the enhanced constant G . To lowest order 
in p the parameter t is also the proper time of comoving matter in its own 
metric since the two metrics coincide under these circumstances. It follows 
that Bubble's constant Hq = A/ A and the current deceleration parameter is 
go = —AA/A"^ and that the critical density for a flat universe is 

So Pc ill our theory is between one tenth and one twentieth the value appro- 
priate to standard General Relativity. This brings it within range of visible 



matter pOf. In fact recent measurements of distant supernovae lead to 



an estimate of go that suggests p = 0.3 x pc and pcc = 0.7 x pc (assuming 



a fiat universe with the further implication that in the present theory, 
visible matter may support the current expansion of the universe albeit with 
the help of a cosmological constant. 



11 Conclusions 

In this paper we have constructed a modified theory of gravity that fits all the 
standard tests of General Relativity that can be made on the scale of the solar 
system including the bending of light by the sun, time delay measurements 
and the precession of Mercury's orbit. The theory is a bi-metric theory of a 
novel kind with a very geometrical structure. It has the flexibility to permit 
the introduction of a galactic length-scale of roughly 30 kpc. Gravitational 
effects of matter at distances below the galactic scale, in the solar system 
for example, are of a conventional kind and of a strength determined by 
Newton's constant, G^ . This outcome is achieved in the theory as the 
result of a competition between two gravitational effects, a strong attraction 
and a repulsion that is nearly as strong. Over the range of the galactic 
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scale the theory allows the repulsion to fade out exponentially leaving the 
much stronger underlying gravity to show through. The gravitational effects 
of matter on this large scale are also conventional in character but of a 
strength determined by the much stronger underlying gravitational constant 
G~ 10 — 20xG^ . These results are consistent with observations of galactic 
clusters and also apply to the expansion of the universe as a whole. As 
a result the critical density for a flat universe is between one tenth and 
one twentieth of the value calculated in General Relativity. This makes it 
more plausible that visible baryonic matter can support the expansion of 
the universe without dark matter, though supplemented by a cosmological 
constant. 

The theory provides a natural explanation for the tendency of dark matter 
to follow the distribution of visible matter since the former is simply an 
amplification of the latter. Dark matter does not have a separate dynamics 
of its own. On this basis the apparent dark matter is not in any way biased 
in its fluctuation structure relative to visible matter. 

We showed that at least for the galaxy NGC3198 the theory can repro- 
duce the rotation curve out to 30 kpc (the end of the measured range) equally 
as well as dark matter models. The theory also predicts that the rotation 
curve will fall away in the 100-200 kpc range and eventually tend to a stan- 
dard inverse square root of distance behaviour but with an apparent mass 
10 - 20 times the expected mass of the galaxy (~ 3 x W^Mq for NGC3198) 
if the conventional Newtonian constant is used for the estimate. Our theory 
however assigns the expected mass but uses the enhanced value for the grav- 
itational constant appropriate at extra-galactic distances. This phenomenon 
should appear for all galaxies. It is therefore important for the theory to test 
rotation curves at distances of 100 — 200 kpc from the galactic centre. It 
is also important to test the effectiveness of the theory in a range of galax- 
ies and to arrive at a consistent picture with a common set of gravitational 
parameters. This task requires a good understanding of the structure of in- 
dividual galaxies and and accurate knowledge of their distance. However we 
note that the establishment of a fixed scale associated with all galaxies would 
be of considerable help in establishing galactic distances. 

In addition to the incomplete phenomenological analyses discussed above 
there are many issues yet to be explored in our bi-metric theory. Of par- 
ticular importance are the development of density fluctations, dynamics and 
formation of galaxies using the modified gravitational law, and evolution of 
the early universe. In this last context the ability of bi-metric theories to 
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support anomalous propagation of signals (VSL phenomena) will be of great 
importance. The existence within the theory of black holes is also a topic 
that it would be very interesting to resolve. On intuitive grounds it seems 
reasonable to suppose that black holes will exist in our theory but the nature 
of the horizon may be more complicated than in General Relativity because 
of its double light-cone structure. 
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Rotation Curve for NGC3198 
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Figure 1: Data for NGC3198 compared to exponential disk model with enhanced 
gravity (full line) and the Newtonian (no enhancement) curve for the same mass 
(dashed line). 
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Extended Rotation Curve for NGC3198 
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Figure 2: The extended rotation curve for the exponential disk model with en- 
hanced gravity (full line) and the Newtonian (no enhancement) curve for the same 
mass (dashed line) 
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